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Abstract 

We examine the violation of the Pauli exclusion principle in the Quasipar- 
ticle Random Phase Approximation (QRPA) calculation of the two-neutrino 
double beta decay matrix elements, which has its origin in the quasi-boson 
approximation. For that purpose we propose a new renormalized QRPA with 
proton- neutron pairing method (full-RQRPA) for nuclear structure studies, 
which includes ground state correlation beyond the QRPA. This is achieved 
by using of renormalized quasi-boson approximation, in which the Pauli exclu- 
sion principle is taken into account more carefully. The full-RQRPA has been 
applied to two-neutrino double beta decay of Ge, 82 5e, 128 Te and 130 Te. 
The nuclear matrix elements have been found significantly less sensitive to the 
increasing strength of particle-particle interaction in the physically interesting 
region in comparison with QRPA results. The strong differences between the 
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results of both methods indicate that the Pauli exclusion principle plays an 

important role in the evaluation of the double beta decay. The inclusion of 

the Pauli principle removes the difficulties with the strong dependence on the 

particle-particle strength g pp in the QRPA on the two-neutrino double beta 

decay. 
23.40.Hc 
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I. INTRODUCTION 



Neutrinoless double beta decay (0u(3j3), which involves the emission of two electrons and 
no neutrinos, demands neutrino to be a Majorana particle and have a non-zero mass. This 
process violates the lepton number conservation and occurs in some theories beyond the 
standard model ( JTJ- for reviews). But, the 0^/9/9-decay has not yet been observed. The 
experimental lower limits on the half-lives provide the most stringent limits on the effective 
mass of the neutrinos and the parameters of the right-handed currents, depending on Ovftft- 
decay nuclear matrix elements. Two neutrino double beta decay (2u/3f3), which involves 
the emission of two antineutrinos and two electrons, is a second order weak process fully 
consistent with the standard model Jl|- ||. It is the rarest process observed so far in the 
nature. Direct measurements of 2i//3/3-decay half-lives provide a sensitive test of nuclear 
structure models. The reliable evaluation of the 2i//?/?-decay nuclear matrix elements is 
necessary to gain confidence in the calculated 0^/3/?-decay nuclear matrix elements. 

The quasiparticle random phase approximation (QRPA) is the nuclear many-body 
method most widely used to deal with the nuclear structure aspects of the double beta 
decay process ||- (T^j. The QRPA has been found successful in explaining the quenching 
of the 0^/9/?-decay nuclear matrix elements and bring them into closer agreement with ex- 
perimental values. But despite this success, the QRPA approach to double beta decay has 
some shortcomings. The main problem is that the results are extremely sensitive to the 
renormalization of the particle-particle component of the residual interaction, which is in 
large part responsible for suppressing calculated two neutrino decay rate. This feature has 
been confirmed in the calculations with both schematic contact 5 interaction and realistic 
finite range two-body interaction. Several alterations of the QRPA have been proposed that 
might change that behavior including proton-neutron pairing |T7|[18| , higher order RPA cor- 
rections QT9J and particle number projection However, none of these extension of 
the QRPA have changed the rapid variation of the calculated 2i//3/3-decay nuclear matrix 
elements with the increasing strength of the particle-particle force. 
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The crucial point of the QRPA consists in the so-called quasi-boson approximation 
(QBA) leading to a violation of the Pauli exclusion principle. The validity of the QBA 
in the evaluation of the nuclear many body Green functions for double beta decay transi- 
tions is questioned because of the generation of too much ground state correlations with 
increasing strength of the particle-particle interaction leading to a collapse of the QRPA. 
The renormalized RPA p2| - [P6[| being proposed many years ago is a method going one step 



beyond the QBA. There are many variations of this method |27|]- |3T|], which basically con- 
sists in replacing the uncorrelated ground state (Hartree-Fock ground state or BCS ground 
state) by the correlated ground state ( RPA or QRPA ground state) in the calculation of the 
expectation value of the commutator of two bifermion operators. The topic of this paper is 
to propose a renormalized QRPA with proton-neutron pairing and to study the influence of 
the ground state correlations beyond the QRPA on the value of 2z//3/3-decay matrix elements. 



We note that recently Toivanen and Suhonen have proposed a similar method |]32] to 
study the nuclear double beta decay which however does not take into account proton- 
neutron pairing. Their method has been applied to calculate the two-neutrino double beta 
decay of 100 Mo by considering some of the ground state correlations beyond the QRPA. 



II. FORMALISM 



We describe the smearing of the Fermi surface by the Hartree- Fock- Bogoliubov 
(HFB) method, which includes proton-proton, neutron-neutron and proton-neutron pair- 



ing 0013 ■ Particle (c+ ama and c rama , r = p,n) and quasiparticle (a+ arria and a^ ama , 



fi = 1,2) creation and annihilation operators for spherical shell-model states labeled by 
(a m a ) are related to each other by the Bogoliubov transformation: 
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where the tilde ~ indicates the time reversed states a Ta m a — (— l) J ' a ~ m ° l a Ta _ ma . The label 
a designates single particle (quasiparticle) quantum numbers (n a ,l a ,j a ). The occupation 
amplitudes u and v and the single quasiparticle energies E aa are obtained by solving the ad- 



equate HFB equation |IT| , |T8"| , |3"3f . The vacuum of the quasiparticle operators will be denoted 
by \0~H FB >. We note that in the ground states of spherical even mass nuclei one can have 
only J=0 and T=l proton- neutron pairs in the framework of the BCS and more general 
HFB theory. In principle one could have T=0, J=odd correlations but only for deformed 
nuclei ||34j| . The treatment here is restricted to spherical nuclei. Such deformation effects 



are taken effectively into account by renormalizing according to [17| , 18| the T=l, J = + 
proton-neutron pairing to be stronger than requested by the Brueckner reaction matrix el- 
ements of the Bonn potential. This strength is fitted to the experimental proton-neutron 



pairing gap extracted from the masses [17,18 



In the limit in which there is no proton-neutron pairing u 2p = v 2p = u ln = v ln = 0. Then 
the HFB transformation in Eq. (|1]) reduces to two BCS transformations, first for protons 
Oip = u P i v i P = V P ) and second for neutrons (u 2n = u n , v 2n = v n ). 

We assume the nuclear motion to be harmonic and the phonon operator to be of the 
form lliTLirglBl 



Q?m* = EK(^ h J n )Af 2 (k, I, J, M) + Y£(k, I, r)A 12 (k, I, J, M)} 



hi 



+ E ^ J *) A U k > h ^ M) + Y™(k } I, J^A^k, I, J, M)h (2) 

k<l 
fj,=l,2 



Here, X m and Y m are forward- and backward- going free variational amplitudes. The 
operator A+ U (k,l, J, M) (/i, v = 1,2) creates a pair of quasiparticles coupled to angular 
momentum J with projection M, namely 

A~^ v {k, /, J, M) = n(kfM, lis) ^2 ^ femfc j;m; a /ijfcm fe a i4m;J 



A^ u (k, I, JM) = (-l) J - M A^(k, I, J, -M) (3) 



The factor n(kfi, Iv) is a normalization in the case in which both quasiparticles are in the 
same shell: 



n(kfi, lv) = (1 + (-1) J 5 kl 5^)/(1 + SuS^f' 2 . 
By using the machinery of the equation of motion we get eigenvalue equation 
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(5) 



with the matrices 

AjOik, vl- y'k', v'l') = (OtpAA^k, /, J, M), [H, AJ V (A/, I', J, M)]]|0+p A ), 

i/Z; //*', v'l') = K PA \{A^(k, I, J, M), [A^tf, 4 M )> H]]\0+ PA ), 
and the unitary-matrix IA 

U = K PA \[A^(k, I, J, M), A+ V ,{k', V, J, M)]|0+, A >. 



(6) 
(7) 



(8) 



Here, if is the quasiparticle hamiltonian and flj^ is the excitation energy — E . and 
Eq are respectively eigenenergies of the excited state | JM^ > (| JM£ > = Qjm*\ q rpa >) 
and of the ground state |0^ PA > {Q^mA^rpa >= 0)- 

In order to solve the Eq. (||) we need additional approximation. In the framework of 
QBA 

{A^(kJ,J,M),A+Ak',l',J,M)} 

* KfbHA^K l,J,M), A i>Ak', V, J, M)]\0+ FB ) 



(9) 



the unitary-matrix U in Eq. @ is just unity-matrix and the matrix equation fl5|) simplifies 



to QRPA with proton-neutron pairing equation (full-QRPA) [[T7| JT8| , |33 1 . 

It is well-known that the QBA violates the Pauli exclusion principle because we have 
neglected terms coming from the commutator. This deficiency could be improved with help 
of renormalized quasi-boson approximation (RQBA) 
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[A^uik-, I, J, M), A+, v ,(k' , I', J, M)] 

^ KpaU^K /, ^ M), A+ U ,(k', I', J, M)\ |0+ PA > 
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If D 



fj,k,vl;J 7 ' 



1, one recovers the usual boson commutation relations for the QBA. The 



reduction of U^ui-j* fr° m 1 reflects the fact that one has the commutation relations of 
Fermion pairs (Pauli principle) and not of bosons. 
It is useful to introduce renormalized operators 

A^(k, I, J, M) = V-^Xj, A+ (*, /, J, M). (12) 

In the RQBA the overlap matrix U is also diagonal: U = D 1 I 2 \ 1 ° \D 1 / 2 . After the 

V° V 

transformations 



X m = V l ' 2 X m , r = V l ' 2 Y m , 
for the phonon operator we can write 

Q3£-= E {x^(kj^nX,(kjA,M)+r^ u (kjAnX u (kj,j,M)} 

fik<i/l 



(13) 



(14) 



and the renormalized quasiparticle random phase approximation (full-RQRPA) equation 
takes the form of RPA-like equation 



(15) 
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The matrix elements A and B are explicitly given by 
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n(kfi, lv)n{kfL, lu){ (E kfl + E^S^S^S^ - (-1)**+* J 5 ki 5^6 lk 6 u/1 ) 

x I G{kal^l5J){u km u lv pu k(ll u iM + v^v^v^v^) 

+F(kal^kjt8J)(u km vi^u^v (i>s + v kfia u lvP v^u (i>s ) 
-(-iyk+^r J F(kal(3(SkjJ)(u km vi^u^ 5 v k ^ + v k ^u^ ] }, (16) 



Bj^fik, vl\ n'k', v'V) = 2 n(kfi, Iv) n(kfa Iv) V][ k 2 >vl . JV V 1 ^, v , v .^ 

x { [ G(kal(3kjlSJ) {u km ui v pv kM v iilS + v klJXX v lv pu kM u iil8 ) 

a/3Sf 

-Fikaipk-ylSJXu^aVivpv^Ufcg + v klwl ui vP u kM v Us ) 
+ (_l)4+ir J ' F(kal(3l5k 1 J)(u km vi vl3 v i , s u kM + v ktU3l ui v pu^ s v k(jn ) ] }. (17) 

Here, G(aabj3cyd5J) and F{aab(3c~^d5J) are respectively particle-particle and particle-hole 
interaction matrix-elements denned in ||35|| . The greek letters a, /3, 7 and 5 are isospin 
indices (e.g., a = p,n) while the greek letters /1, //, ^ and v' run the value 1 or 2 and denote 
the "quasi-isospin" of the quasiparticles. 

In order to calculate the RQRPA matrices Aj-* and B J7 , one has to determine quantity 
T^ink^i-^ ■ We follow the method of ref. |3T| and express the one-body densities of Eq. ( |lCf ) 



in terms of mappings 

Wb = i E z ' J ' M ) J, 4 m). (is) 

By using eqs. (HD, fllj), (U), (UD together with 

XviK 1, J, m) = J2(x;,(k, 1, r)Qj&, - F£(*, /, J*)G3W, (19) 

m 

we obtain the following expression for V^^i-j^: 
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DkulvJ^ -1-^ E ^fc^VJ'^'^I^M^'^' k'> ^ )| 2 



* E ^v^lOU',^')! 2 (20) 



Ji 



I'w 1 



The system of non-linear equations in Eq. (P0|) could be solved numerically by an iteration 
process. As input the solution of the full-RQRPA equations for all the considered multi- 



polarities is required. However, in order to calculate the full-RQRPA equation in Eq. (JT5 
we need the knowledge of the renormalization factors T>k^i v j^- Therefore, the selfconsistent 
scheme of the calculation is a doubly iterative problem. As starting values for T> k ^ v j^ for 
the first iteration could be used e.g., 

ZW = l-4 E j' 2 \Y™,(k,k',J'«')\ 2 -± £ j' 2 \Y™,(l,l',J'*')\ 2 , (21) 

Jk fcV h I 1 " 1 

with the amplitudes F m being the solution of the full-QRPA equation. 

The half-life for 2z//3/3-decay could be written in a factorized form |T[- |J: 

(T'yfy 1 = G% T \M% T \\ (22) 

where G 2 q T is a lepton phase space integral and the nuclear matrix elements M GT in the 
full-RQRPA is 

_ / < 0^ PA ||r+a||l+ >< 1+ |1+ >< 1+.\\t + v\Kpa >i , s 

M 2 G t= E " f S ' . ' (23) 



i' f 



with reduced transition-amplitudes of the Gamow- Teller beta decay operator 
<l^ii-r + tr||0+ PA > i = E *(M) PUi+) 1/2 

/ik<ul 

x [ X^ik, I, l + )(u l kfip v l lvn — v l k ^ n u\ up ) + Y^(k, 1, l + )(v k/Jip u % lvn — u % k ^ n v\ up ) ], (24) 



/ U/2 
k/ilvl+J 

X [ -^"iw/ ^3 l + )(' y fc M p' lt Ln — U kfm V Lp) + Y ^{k, I, ^- Jr ){u I klxp Vi vn — V klln u{ vp ) }. (25) 



/<0^||r + a||l+ / > = ^ <KM)(2>Ui + 



Here, a(k,l) =< k\\a\\l > is the reduced matrix element of a operator. The HFB u and 
v factors entering in eqs. fl24 ) and (|25|) belong to two different even-even system. The 



occupation amplitudes u l ,v l and u^,v^ are obtained by solving the HFB equation for the 
initial and final nuclear states, respectively. |1+ > and |1+ > are respectively the excited 
states calculated from the ground state of initial (A,Z) and final (A,Z+2) nuclei with the 
excitation energies and The states |1+ > (|l£ >) are expressed in terms of 

forwards- ~X%(k,l, 1+) (X%(k, 1, 1+)) and backwards- Fj(M,l + ) (F™f (jfc,Z, 1+)) going 
amplitudes. Unfortunately the states |mj, 1 + > and \rrif, 1 + > are not orthogonal to each 
other. For the overlap matrix of these states we write 

We note that in the limit of small ground state correlations T>k tl wj^ — 1 the full-RQRPA 
could be replaced by the full-QRPA. The deviation between the results of both methods is 
a signal showing on the restriction of the validity of the QBA in a given physical region. 
Further, in the limit in which there is no proton-neutron paring, the full-QRPA reduces 
to two matrix equations The first pp+nn QRPA equation has the origin in 



the proton-proton and neutron-neutron quasiparticle excitations. The second pn-QRPA 
equation considers only the proton-neutron quasiparticle excitations from the ground state. 
By neglecting the proton-neutron pairing interaction the full-RQRPA equation decomposes 
in a similar way on pp+nn RQRPA and pn-RQRPA equations. However, pp+nn RQRPA 
and pn-RQRPA equations remain coupled through the Eq. ( p0|) . Therefore, both equations 
have to be solved parallel in the iterative process. 



III. CALCULATION AND DISCUSSION OF THE RESULTS 

We applied the full-RQRPA method to the 2z//5/3-decay of 76 Ge, 82 Se, 128 Te and 130 Te. 
We assumed the single particle model space both for protons and neutrons as follows, 
(i) For 76 Ge — > 76 Se and 82 Se -^ S2 Kr the model space comprises 13 levels: 

1^1/2; 0^5/2) 0^3/2) 1^3/25 1^1/2) O/7/2 > 0/5/2 5 ^ S 1/2j 1^5/2) 1^3/2) %9/2) 0fl , 7/2) 0^11/2- 

(i) For 128 Te ^ 128 Xe and 130 Te ^ 130 Xe we used 16 levels: 
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1^3/2) lPl/2) 0/7/25 0/5/2! 2s 1//2 , ld 5 /2, 1^3/2) 0^9/2, 0g 7 / 2 , 0/i 11//2 , 0/l 9 / 2 , 1/7/2) 1/5/2) 



2^3/2) 2px/ 2 , 0i 13 / 2 . 

The single particle energies have been calculated with a Coulomb-corrected Wood-Saxon 
potential. For the two body interaction we used the nuclear G-matrix calculated from Bonn 
one-boson exchange potential. The single quasiparticle energies and occupation amplitudes 
have been found by solving the HFB equation with p-n pairing for both the parent and 
the daughter nuclei in the above mentioned space. The renormalization of the proton- 
proton, neutron-neutron and proton-neutron pairing interaction has been determined ac- 



cording to ref. pj|,^o| . The renormalization of the T=l J=0 proton-neutron force is due to 
the quadrupole deformation of the intrinsic nucleus. It could also be described in a spher- 
ical basis by admixing to the single particle propagator 2 + quadrupole excitations. Thus 
this renormalization and the HFB step influence only the quasiparticle states. The collec- 
tive excitations are described in a second QRPA step. The renormalization we fit to the 



proton-neutron pairing gap ]T^,[TS[ . 

In the calculation of the full-RQRPA equation we renormalized particle- particle and 
particle- hole channels of the G-matrix interaction by introducing two parameters g pp and 
g ph , which in principle should be equal to unity. Since the matrix element Mq^ is much 
more sensitive to g pp than to g ph , we will set g ph = 1.0 and discuss the dependence on 
g pp . Only for the 2 + -channel of the interaction g ph was fixed to g ph = 0.8 for all studied 
2i//3/3-decay transitions. For higher value of g ph the particle-hole interaction in 2 + channel 
is too strong. The lowest eigenvalue becomes the imaginary and leads to a collapse of the 
correlated ground state. 

It is worthwhile mentioning that the calculation of Mq u t within full-RQRPA needs a 
great computational effort. We remind that the full-RQRPA self- consistent scheme of the 
calculation requires the solution of the full-RQRPA equation (|T5| ) for all multipolarities J w 
in each iteration for the initial and the final nuclei. In comparison with the full-RQRPA, 
the full-QRPA calculation of M GT requires to solve the full-QRPA equation only for the 
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multipolarity 1 + once for initial and once for final nucleus. The iterative procedure of the 
full-RQRPA have been found to converge rapidly. 

As the full-RQRPA calculation is a very time consuming process, it is interesting to 
determine, which multipolarities J n are important to take into account to calculate the 
renormalization factors V^^ij* in Eq. fl2"0"|). We made a calculation for g pp = 1.0 including 



more and more multipolarities J n in the renormalization scheme of the full-RQRPA. The 
results are shown in Fig. 1. Reading from the left to the right, the quoted multipolarity is 
added at every point to the ones of the preceding point. One can clearly see a saturation effect 
for large multipolarities. This can be explained by the fact that the higher multipolarities 
are less collective and do only little contribute to the ground state correlations. Practically, 
one is able to find a combination of the multipolarities, that optimizes both the effect of 
renormalization and the effort of the computing-time. Fortunately, it turns out that one 
combination, which consists of the 1 + , 2 + , 2~ and 3~ multipolarities, shows as good results 
as the whole set. The corresponding values of M"q T for the studied nuclear transitions are 
plotted with empty symbols in Fig. 1. A given combination of the multipolarities have been 
used to study the g pp dependence of the matrix element Mq T . 

The goal of the present paper is to study the reliability of the QBA in the nuclear 
structure calculation of M"q T . This information can be obtained by the comparison of the 
results of the full-RQRPA and the full-QRPA in the physically acceptable region of the 
g pp parameters 0.8 < g pp < 1.2. Our results for 2v(3(3- decay of 76 Ge, 82 Se, 128 Te and 
130 Te are presented in Fig. 2, Fig. 3, Fig. 4 and Fig. 5, respectively. We see that the 
effect due to the ground state correlations beyond the RPA is improving the agreement 
for all studied nuclear transitions. The inclusion of the ground state correlations beyond 
RPA was found to stabilize the behaviour of M^Y T as a function of g pp . The full-QRPA 
and the full-RQRPA equations are of the same form, differing only by a renormalization 
of the two-body interaction by the factor V^k,ui-,j^ [see Eq. (p~5|)1. For V^k,ui-,j^ — 1, what 
is the case for small ground state correlations, the results of both methods are the same. 
However, the QRPA quenching mechanism for M(Y T has its origin in generating too much 
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ground state correlations with increasing strength of the particle-particle interaction. Strong 
ground state correlations lead to a significantly smaller value of V ^k,vi-,j-" and reduces the 
strength of two-body interaction. In this way a damping of the dependence of Mq u t on g pp is 
obtained. The full-RQRPA has been found to reproduce well the value of M^f exp deduced 
from the experimental half-live. Agreement is achieved for w Ge and 82 Se for the parameters 
g ~ 0.95 and g pp ~ 1.00, respectively. In the case of 2i//3/3-decay of 128 Te and 130 Te we 
need higher values of the renormalization parameter g pp : g pp ~ 1.3. We remind that the 
Te isotopes are known to be a very special and delicate case. W. Haxton already found 



this out [|36]]. It could be connected with the fact that the two Te isotopes 128 and 130, for 
which the double beta decay is measured, are unharmonic or even slightly deformed, so that 
a spherical RPA treatment as used here is not sufficient. 

Both, the full-QRPA and full-RQRPA, include proton-neutron pairing correlations via 
the HFB method. In order to clarify the role of proton-neutron pairing in the stabilization 
of the dependence of Mq u t on g pp , we switched off proton-neutron pairing interaction and 
performed the calculation within pn-QRPA and pn-RQRPA (no p-n pairing). The results are 
presented in Figs. 2-5. We can see that also in this case one can avoid a strong instability of 
Mq u t . It is interesting to note that the effect due to the proton-neutron pairing is significant. 
The inclusion of proton-neutron pairing in the self-consistent scheme of the calculation is 
also helping to shift the instability of M"q T on g pp to larger values of g pp , which are outside 
the physical range (0.8 < g pp < 1.3). 



IV. CONCLUSIONS AND OUTLOOK 



In summary, we have proposed a new higher (renormalized) QRPA many-body method 
for nuclear structure studies which we call full-RQRPA. We point out that in the full-RQRPA 
proton-neutron pairing is taken into account. The full-RQRPA includes the Pauli exclusion 
principle in a more proper way in comparison with the QRPA by using the renormalized 
quasi-boson approximation. The one-particle densities of the correlated state have been 
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evaluated by following the approach of Catara et al. In this way the ground state correlations 
beyond the RPA are included in the self-consistent scheme of the calculation consisting of a 
set of non-linear equations in amplitudes X and Y, which are solved by iteration. 

We have employed the full-RQRPA to calculate 2i//3/3-decay of 76 Ge, 82 Se, 128 Te and 
130 Te. The nuclear matrix element Mq u t governing this process has been calculated. We 
have found M"q t to be rather stable in the physically acceptable region of the parameter 
g , which determine the strength of the particle-particle interaction. It allows us to predict 
more reliable values of the nuclear matrix elements. The relative large differences between 
the results of the full-RQRPA and the full-QRPA are related to violations of the Pauli 
exclusion principle by the second method due to ground state correlations. The quasi- 
boson approximation is a poor approximation to study the 2i//3/?-decay. Therefore, an 
alternative methods should be elaborated in which the Pauli exclusion principle is explicitly 
incorporated. The shell model approach, which does not violate the Pauli exclusion principle, 
fails to construct all needed states of the intermediate nucleus. We note that the recently 
proposed two-vacua random phase approximation method could be perhaps a step in the 
right direction pT7| , |38[ . 

The renormalized QRPA which includes the Pauli effect of fermion pairs goes beyond the 
quasi-boson approximation. The inclusion of the Pauli principle reduces the ground state 
correlations and moves the collapse outside the physical region. Our matrix elements are 
derived from the Bonn-potential and if we would stick to them, g pp would be unity. Each 
RPA solution is collapsing eventually if the forces are made stronger and stronger. The 
point is that this collapse happens here only outside the physical region as in all stable RPA 
solutions and therefore the collapse disappeared for the description of the nuclei. 

The full-RQRPA could find a wider use in the studies of other nuclear processes, e.g., 
single beta decay, neutrinoless double beta decay and pion double charge exchange reactions. 
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FIGURES 
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FIG. 1. The influence of including more and more multipolarities J 11 in the full-RQRPA 
calculation of Mq T for g pp = 1.0. Starting from 1 + the mentioned multipole is added to the former 
ones at every point. By this procedure the last point shows the value of M"q T obtained by the 
consideration of 12 multipolarities ranging from + to 5 _ . The empty symbols mark the values for 
the combination of 1 + , 2 + , 2~ and 3~ multipolarities, which turn out to exhaust approximately 
the most of the correlations beyond RPA. 
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76 Ge -> 76 Se 
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- full-RQRPA 




'PP 



FIG. 2. The Gamow- Teller transition matrix element Mq^ of the 2z//3/3-decay of Ge is plot- 
ted as function of particle-particle coupling constant g pp . The solid line corresponds to full-RQRPA 
(with p-n pairing), the dashed line to full-QRPA (with p-n pairing), the dash-dotted line to 
pn-RQRPA (without p-n pairing) and the dotted line to pn-QRPA calculation (without p-n pair- 
ing), respectively. 
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82 Se -> 82 Kr 
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FIG. 3. The same as Fig. 2 for 82 Se. 
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FIG. 4. The same as Fig. 2 for 128 Te. 
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FIG. 5. The same as Fig. 2 for idU Te 
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